Light Deflection, Lensing, and Time Delays from Gravitational Potentials and 
Fermat's Principle in the Presence of a Cosmological Constant 
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The contribution of the cosmological constant to the deflection angle and the time delays are 
derived from the integration of the gravitational potential as well as from Fermat's Principle. The 
findings are in agreement with recent results using exact solutions to Einstein's equations and 
reproduce precisely the new A-term in the bending angle and the lens equation. The consequences 
on time delay expressions are explored. While it is known that A contributes to the gravitational 
time delay, it is shown here that a new A-term appears in the geometrical time delay as well. 
Although these newly derived terms are perhaps small for current observations, they do not cancel 
out as previously claimed. Moreover, as shown before, at galaxy cluster scale, the A contribution 
can be larger than the second-order term in the Einstein deflection angle for several cluster lens 
systems. 

PACS numbers; 95.30.Sf,98.80.Es,98.62.Sb 



I. INTRODUCTION 

Various complementary observations seem to indicate 
that the expansion of the universe has entered a phase 
of acceleration 1] . This has become one of the most im- 
portant current problems in cosmology and all physics, 
as discussed, for example, in and references therein. 
While several explanations are possible current ob- 
servations 1] are consistent with a cosmological constant 
in the Einstein Field Equations. 

In a recent work, the authors of 0] used the 
Schwarzschild-de Sitter (SdS) metric and showed that 
the cosmological constant does contribute to the bend- 
ing of light around a concentrated mass. It was shown 
there that even though A drops out from the null geodesic 
equation, it still contributes to the deflection of light sim- 
ply because of the geometry of the SdS spacetime repre- 
sented by the metric, regardless of other considerations. 
The result made a correction to the ques tion (see e.g. 
@,i,0,[iji) and was confirmed in Next, in 

reference [J| , the authors brought the result into an obser- 
vational context and used an exact solution construction 
where a Schwarzschild-de Sitter vacuole was embedded 
into a Friedmann-Lemaitre-Robertson- Walker (FLRW) 
background, and where a A-term was added to the deflec- 
tion angle within the broadly used lens equation. Finally, 
using observations of Einstein angles around clusters, the 
authors were able to put upper-bounds on the value of 
the cosmological constant, only two orders of magnitude 
away from the value determined by cosmological probe 
constraints. 

In this paper, we use the integration of the gravita- 
tional potential method and also Fermat's Principle to 
calculate the contribution of the cosmological constant 
to the deflection of light, the lens equation, and the time 
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delays. Whereas it is known that A contributes to the 
gravitational time delay, we show here that a new A-term 
appears in the geometrical time delay. 

First, It is worth clarifying a point that has been a 
source of ambiguity in some of the recent literature about 
the effect of A on the bending of light and gravitational 
lensing by a galaxy cluster in an FLRW spacetime. As we 
show in Figure 1, the light bending due the lens occurs in 
a region close to the lens (SdS vacuole) and then once the 
light transitions out of the vacuole into FLRW spacetime 
then the bending stops. Away from the cluster (i.e. out- 
side the vacuole) , the ray of light propagates in an FLRW 
background. Now, from the geometry, the lens equation 
includes the bending angle and also angular diameter dis- 
tances. While A is present in the lens equation via the 
usual expressions for the angular diameter distances in 
FLRW, A also enters the lens equation via the bending 
angle itself that occurs inside the SdS vacuole. As evalu- 
ated in pj , the contribution of A to the deflection angle is 
small but larger than the second-order term in the Ein- 
stein angle for several cluster lens systems. Moreover, 
this is also a matter of mathematical correctness since 
the newly derived terms do not cancel out as previously 
claimed. In relation to this clarification about the lens 
geometry above, recent papers [ll|, [l^l refer to our pre- 
vious work [j] without taking into consideration the 
construction described above and used in [5| . Then, it is 
concluded in ^13i] that terms of the order of ^ ArgVQ (in 
their notation) may exist. That term is precisely the A- 
term derived in our previous work [^]. It is hopped that 
this paragraph along with Figure 1 will clarify this point 
about the cosmological constant and the lens equation 
geometry. 

It was recently pointed out in that some of the A- 
terms contributing to the deflection angle could be incor- 
porated into the angular diameter distances when these 
are extended up to the lens and not the boundary of the 
vacuole while other A-terms cannot; this will be investi- 
gated elsewhere. 
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FIG. 1: The lens equation geometry. Observer, lens, and source are at O, L, and S, respectively. The position of the unlensed 
source is at an angle /3, the apparent position is at the angle 9 and the deflection angle is a. The distance from the observer to 
the source is Dos, from the observer to the lens is Dql, and from the lens to the source is Dls- As usual, the lens equation 
follows from the geometry as 6 Dos = PDos + oiDls- A ray of light propagates from the left within the FLRW background, 
then enters inside the Schwarzschild-de Sitter vacuole around the lens where it gets deflected by the lens, and then it exits the 
vacuole and continues its propagation in the FLRW background to the right. The results derived from matching exact solutions 
to Einstein's equations in [j| and those derived in this work from the usual approximate constructions used in gravitational 
lensing literature are in perfect agreement for the contribution of A to the deflection angle and the lens equation. The same 
result is also derived using Fermat's Principle. 



It is also worth noting that the contribution of the cos- 
mological constant to the bending angle enters via the 
Weyl focusing of the Sachs equations [H, [H, [O] and 
takes place inside the SdS vacuole (see also footnote [22] 



ground and was given by 
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while outside the vacuole only the 
id \vf\ takes place and the cosmo- 
logical constant appears via the angular distances in the 
FLRW background. This applies to the bending angle 
and the related time-delay as well. 

As we will show, the A-term in the lens equation as 
discussed in [3| can be derived using various methods 
used in gravitational lensing and propagates to time delay 
calculations. 



II. DEFLECTION POTENTIALS, THE LENS 
EQUATION, AND THE COSMOLOGICAL 
CONSTANT 



In this section, we show that the contribution of 
the cosmological constant to the light bending angle 
and the lens equation \a] can also be derived from 
gravitational potentials, a method that is frequently 
used in gravitational lensing literature, see for example 

[illil[il,|2i. 



Ishak, Rindler et al. derived in [4| an amended lens 
equation where the contribution of the cosmological con- 
stant was added to the Einstein bending angle. The 
A-term was derived there using an exact construction 
of a Schwarzschild-de Sitter (SdS) vacuole embedded in 
a Friedmann-Lemaitre-Robertson- Walker (FLRW) back- 
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(1) 



where i? is a radius close to the lens as used in [3, 0| 
and similar to the impact parameter, e.g. [13, H^, but 
distinct from it since the SdS spacetime is not asymptoti- 
cally flat, and rf, is the radial coordinate at the boundary 
of the vacuole where the spacetime transitions from the 
SdS spacetime to an FLRW background. 

We note that one can use again the following matc hing 
conditions at the boundary of the vacuole, see e.g. |24 . 
MM, 



n in SdS = a,{t) n in 



FLRW 



and 



47r 



msds 



in SdS ^ Pmatter in F LRW , 



and re- write the Ishak- Rindler A-term, ([T]), as 



"A 



All / Smsds 

3 Pmatter in FLRW 



(1/3) 



(2) 



(3) 



(4) 



This expression gives the A-term as a function of the 
mass of the galaxy or cluster lens in the SdS vacuole and 
the matter density of the FLRW background. 

The second method that we shall employ to calculate 
a\ is analogous to the method of approximation that is 
frequently used in gravitational lensing literature where 
the lens (inhomogeneity) in an FLRW background is rep- 
resented by a Newtonian potential inserted in a post- 
Minkowskian line element or a post-FLRW line element 
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(see for example [1^, [l^, [13, 1201 ) • The metric in such a 
construction is then given by gab — Vab + ^ah where hat 
measure the departure from the Minkowskian metric rjab- 
First, we recah that hat in hnearized Schwarzschild 
spacetime and the associated Newtonian potential $ can 
be read off the usual Schwarzschild metric in isotropic 
coordinates given by (e.g. [ll], |2^) (we use relativistic 
units where we set G — c = 1) 



1 



2r 



1 + — 



m \ 



dt' + {l + —j (dr'+r'dn'), (5) 



and are given by 



htt 



m 
r 



(6) 



For a ray of light traveling in the x-direction, the first- 
order term of the Einstein deflection angle is then given 
by 



V_L(/itt + hxx)dx 



(7) 



where = V — V|| is the gradient transverse to the 
path, see for example [1^, . At the vacuole boundary, 
Xb = yjrf^ — B? and a straightforward integration yields, 



the Einstein first-order term 



4m 



(8) 



For the contribution of A to the bending angle which 
happens close to the lens (i.e. inside the vacuole re- 
gion), some subtleties come into consideration as we shall 
see. In linearized GR, we can add the hab due to vari- 
ous sources. It therefore follows that the bending an- 
gle itself can be obtained by adding the bending angles 
due to different hab- Moreover, in order to comply with 
the isotropic FLRW background outside the vacuole, we 
should here also use isotropic coordinates for the de Sitter 
metric inside the vacuole, see e.g. [13, and reference 
therein. 



ds^ 



1 - 



1 



12 



dt^ 



1 



{dr^+r^dn'^) (9) 



which linearizes to gab = Vab + hab with, see for example 



htt = -2$ = 
hoi = 



Ar2 



Ar2 



hij = -2'I'(5y = — 6"*^*^' 



(10) 



The A contribution to the bending angle is then found 
by using again 



ttA = 



^T. \ V ±{htt + hj;x)dx 
V_L{<i> + ^)dx ■ 



ARn 



(11) 



Now, combining equations ([S]) and (fTT|) , using the lens 
geometry (Figure 1), and the usual assumption of small 
angles, the lens equation reads 



f3 = [aEinstein + Q^a) 



Dls 
Dos' 



(12) 



The A-term ((TT|) and the lens equation ([T^ are exactly 
the A-term and the lens equation as derived in 4] from a 
different method based on the exact matching of an SdS 
spacetime vacuole to an FLRW spacetime background 



In section (IV) below we shall give yet another alter- 
native for deriving the A-term, this time from Fermat's 
Principle. 



III. TIME DELAYS AND THE 
COSMOLOGICAL CONSTANT 

A ray of light traveling in the gravitational field is 
known to experience two types of time delay, see for ex- 
ample [isl [so] ■ The first is the geometrical time delay. 



' geom , 



due to the extra path length resulting from the 
deflection. The second is the gravitational time delay, 
Tgrav, due to the grav itational potential, known as the 
Shapiro delay [isl l29l|. The contribution of the cosmo- 
logical constant to the gravitational time delay has been 
discussed in for example [3, [13- However, we show here 
that A also affects the geometrical time delay because of 
its contribution to the bending angle shortens the bent 
path. 

Before we derive the two terms, we want to clarify that 
the contributions of A that we are interested in here are 
distinct from the obvious ones that come via the angular 
diameter distances as we explained in the introduction, 
but we are interested in the contributions that happen 
inside the vacuole around the lens. 

First, the gravitational time delay due to A can be 
derived using (jTUj) , see for example [27[ , and is given by 



TA,, 



Ar| 
18 



{htt + h.xx)dx 
Arb 



12 



(13) 



This result gives explicitly the two leading A-terms for 
the contribution of A to the gravitational time delay. One 
can see that these terms are in agreement with [27| when 
their integration boundaries pi and p2 are set equal to 
+rb = {+Xb,R,Q) and — rb = (— 0), respectively. 
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Now, we recall the Shapiro time delay due to the 
Schwarzschild potential (O. It is given for our vacuole 

by 



'^Shapiro 



^mdx — 4m log 



R 



(14) 



So the gravitational time delay due to A adds to the 
Shapiro delay. These time delays are due to the gravita- 
tional potentials which cause clocks to slow down |2ll.[23j. 

For distant galaxies and clusters of galaxies, the 
Shapiro delay and the A-term need to take into account 
the redshift of the lens, the geometry as given by Figure 1, 
and the angular diameter distances as given in the FLRW 
metric away from the potential. Also, when dealing with 
measurements and observations, the gravitational time 
delays req uire modeling of the lens systems, see for ex- 
ample [131 HH. Following the standard approach [Tsllsot. 
we put the two time delays in a cosmological setting as 

"^grav — (1- ~^ ) r-j ij'Shapiro ^ "^A^grav^ (-^^) 

^ LS 

where all the quantities are as defined in Figure 1 and 
the angular-diameter distance is given by 



D{z) = 



dz' 



(16) 



in a spatially flat cosmology. 

Next, we derive the geometrical time delay, T\^geom^ 
due to A. This term is new and results from the contri- 
bution of A to the deflection angle and thus to the path 
length. A general derivation of Tgeom from the lens geom- 
etry (e.g. figure 1) and the FLRW metric can be found 
in for example |l5|. The geometrical time delay is given 

by 



Pf. (17) 



We can use the lens equation (fT^ to rewrite (fT7|l as 



'^geom — (-^ ~l~ -^l ) 



Dos 2 



(18) 



Now, we can expand to write 



Tgeom ~ (1 + ^t) p; 77 ['^Einstein + (^Einstein OlA) 



where we define the last term as 



(19) 



TA,, 



OiEinstein Q^A 



n J- ^ -Pj-s Dql ( 4A mr,, 
^ D 



(20) 



The term (j20p is new and represents the contribution of 
the cosmological constant to the geometrical time delay 



that enters from its contribution to the deflection angle. 
We clarify again that the contribution of A as it appears 
in the last fraction in (|20p is separate from the contri- 
butions that A has via the angular diameter distances 
as given in the FLRW background, outside the vacuole. 
This newly derived term is perhaps observationally small 
but does not cancel out to zero as previously claimed. 

The sign of Tx.geom is negative. So while the effect 
of a Schwarzschild mass is to produce a retardation due 
to the deflection of the light that it creates, the effect 
of A on the geometrical time delay is the opposite. A 
positive cosmological constant diminishes the deflection 
angle and the length of the path that light needs to travel 
and hence it diminishes the geometrical time delay. 

Finally, we note that when using observations, it is 
rather the relative time delays between different images 
that are used [13, [IH , i.e. 



Ari: 



(21) 



Also, while observations fairly directly yield measure- 
ments of redshifts of lenses and deflection angles, time 
delay measurements require some sophisticated model- 
ing of lens systems [13, [HI ■ These are out of the scope 
of this theoretical work but need to be studied in future 
work using modeling and simulations of observed cluster 
lens systems and images. 

We will show in the next section that the contribution 
of A to the deflection angle and consequently to the ge- 
ometrical time delay can also be derived from Fermat's 
Principle while the Shapiro delay follows directly from 
the spacetime metric. 



IV. FERMAT'S PRINCIPLE AND THE 
CONTRIBUTION OF A TO THE LENS 
EQUATION AND THE TIME DELAYS 

The expression for the deflection angle also follows 
from using Fermat's Principle and the Euler-Lagrange 
equations of the variational Principle 



5 I ndl 



(22) 



where n is considered as an effective index of refraction 
of the gravitational held and dl is the path of the ray of 
light. 

First, as discussed in for example [H, HH, the null 
curve ds^ = gabdx°'dx^ = for the Schwarzschild poten- 
tial ([6]), leads to the Fermat's Principle ((22)l . where 



n = 1 - 2$, (23) 
from which the Euler-Lagrange equations give [l^, [ll| , 

4m 



R 



(24) 



Now, we apply Fermat's Principle in order to derive the 
contribution of the cosmological constant to the bending 
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angle and use hab as given by (fTO|) . For a future- directed 
null curve, 



ds^ = 



'dt^ (1 - 2$) + (1 - 2*) dp, 

+ (25) 



or simply 



dt = (1 - ($ + *))<i? (1 + (26) 

First, the gravitational time delay follows from the sec- 
ond term of the RHS of equation The integral of 
which gives 



At^dl- J{^ + 'i')dl = dl + TA^gr 



(27) 



where TA^grav was already calculated in (jl3p of the pre- 
vious section. 

Next, it follows from equations ([H]) and (|26p that 



n = l-($ + *) = l + — . 



(28) 



Following [1^1 , we consider a ray of light traveling along 
the path dl with unit tangent vector e. The deflection 
angle is given by the change in the direction of the null 
ray. From our equations ([H]) and (|28p and the Euler- 
Lagrange equations, it follows that 



de 



and 



(V($ + *) - e (e . V($ + *)) = + *) 



j V_L($ + *)d?, 



(29) 



(30) 



a result already derived in pip using null geodesies and 
the metric. So for a null ray traveling in the x-direction 
with integration boundaries as given in (jllll . this ap- 
proach gives again the A-term as given in 4] and re- 
derived in section II, i.e. 



"A 



time delay calculations. We used here an approximate 
construction that is frequently used in gravitational lens- 
ing literature and where the lens (inhomogeneity) 

in an FLRW is represented by a potential inserted in a 
post-Minkowskian or post-FLRW line elements. Our re- 
sults for the deflection angle and the lens equation are 
in perfect agreement with previous work '3] where 
an SdS spacetime vacuole was exactly embedded in an 
FLRW spacetime background. Then we derived expres- 
sions for the effect of A on the gravitational and geo- 
metrical time delays. Whereas it is known that the cos- 
mological constant contributes to the gravitational time 
delay, we showed here that a new A-term appears in the 
geometrical time delay. 

We also clarified the geometrical figure of interest 
whether one uses a construction based on exact solu- 
tions to Einstein's equations or based on approximate 
constructions. As shown, in Figure 1, the bending of 
light occurs in the region close to the deflector inside the 
SdS vacuole. Outside the vacuole, the ray of light sim- 
ply propagates in an FLRW spacetime. The contribution 
of the cosmological constant to the bending angle enters 
via the Weyl focusing of the Sachs equations [l^ [l^ [13 
and takes place inside the vacuole, while outside the vac- 
uole only the Ricci focusing [l^, [3, [l7| takes place and 
A appears via the angular distances in the FLRW back- 
ground. 

Indeed, the A-terms that are derived here and in Q 
are related to the bending of light which occurs inside the 
SdS vacuole, close to the lens. These contributions are 
separate from the usual presence of A in the expressions 
for angular diameter distances in the FLRW background. 
It was recently discussed in reference JJi] that some of the 
A-terms contributing to the light bending angle could be 
incorporated into the angular diameter distances when 
these are extended up to the lens and not the bound- 
ary of the vacuole while other A-terms cannot, this will 
be investigated elsewhere. Although the newly derived 
A-terms for the lens equation and the time delay are per- 
haps small for current observations, they do not cancel 
out as previously claimed. Moreover, at galaxy cluster 
scale, it was shown in (3] that the contribution of A to 
the bending angle can be larger than the second-order 
term in the Einstein deflection angle for several cluster 
lens systems. 



(31) 



V. CONCLUSION 

Using the integration of gravitational potentials and 
Fermat's Principle, we find that a cosmological constant 
term appears in the deflection angle and propagates to 
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